Enhanced soliton transport in quasi-periodic lattices with short-range aperiodicity 
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We study linear transmission and nonlinear soliton transport through quasi-periodic structures, 
which profiles are described by multiple modulation frequencies. We show that resonant scattering 
at mixed-frequency resonances limits transmission efficiency of localized wave packets, leading to 
radiation and possible trapping of solitons. We obtain an explicit analytical expression for optimal 
quasi-periodic lattice profiles, where additional aperiodic modulations suppress mixed-frequency 
resonances, resulting in dramatic enhancement of soliton mobility. Our results can be applied to 
the design of photonic waveguide structures, and arrays of magnetic micro-traps for atomic Bose- 
Einstein condensates. 
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Periodic structures can provide an efficient control of 
wave transmission and localization, making it possible 
to realize new physical regimes that are not allowed in 
homogeneous systems. Photonic crystals and Bragg grat- 
ings are used to tailor dispersion, diffraction, and emis- 
sion of electromagnetic waves 0, Q and similar concepts 
have been developed and demonstrated for management 
of atomic Bose-Einstein condensates in periodic poten- 
tials 0, . Quasi-periodic photonic structures such as 
optical super-lattices and 2D quasi-crystals can offer even 
more flexibility in designing theproperties of optical 
and atomic matter waves 0, The modulation of 
quasi-periodic structures is defined by multiple incom- 
mensurable spatial frequencies, resulting in the lack of 
translational symmetry which is a feature of conventional 
periodic structures. This gives rise to an important ef- 
fect of stronger field confinement at certain locations due 
to position-dependent coupling between the lattice sites, 
which can be utilized for various applications including 
efficient lasing 

Many novel phenomena in periodic structures origi- 
nate due to the modification of the linear spectrum which 
is separated into a number of transmission bands and 
band-gaps, where waves experience total reflection 
For a one-dimensional lattice with period the gaps are 
centered at the frequencies corresponding to Bragg reso- 
nance condition for the wavevector component across the 
lattice, 2kx = mK, where m = ±1, ±2, ... is the order 
of resonance and K = 27T/d is the wavenumber of lattice 
modulation. In a quasi-periodic structure, primary gaps 
are defined by several dominant frequencies Kj^ however 
there also appear multiple resonances at mixed frequen- 
cies Kj ± Ki, resulting in a complicated spectrum con- 
taining multiple mini-bands and mini-gaps 9]. This can 
have a detrimental effect on propagation of wavepackets 
through the lattice, which can be inhibited ^ or ac- 
companied by strong radiation . It has been predicted 
that nonlinear wave self-action, which is known to ap- 
pear due to atom-atom scattering in BEC or light-matter 



interactions in optical crystals, can suppress Anderson 
localization and allow wave transmission through disor- 
dered systems |^ . We find that such nonlinear transmis- 
sion may be also possible in quasi-periodic lattices in the 
form of gap solitons (see Refs. [HlllElElEllil, and 
references therein), however they still suffer from radia- 
tion losses that increase for stronger localized wave pack- 
ets which spectra overlap mixed-frequency resonances. 
Resolution of these limitations may open new possibili- 
ties for realizing efficient transport of optical and matter 
waves in quasi-periodic structures. 

In this Letter, we suggest a novel approach to the de- 
sign of quasi-periodic structures. We show that mixed- 
frequency resonances can be suppressed through addi- 
tional aperiodic lattice modulations which inhibit coher- 
ent back-scattering at specific frequencies. Aperiodicity 
if often associated with disorder, leading to the appear- 
ance of critically and fully localized states or Anderson 
localization j^, [l^. However, we present an analytical 
procedure which can be used to define, in a systematic 
way, quasi-periodic structures with short-range determin- 
istic aperiodicity possessing regular (absolutely continu- 
ous) band-gap spectrum, where all eigenmodes are ei- 
ther extended or exponentially decaying. These 'optimal' 
structures allow for unrestricted motion across the lattice 
of solitons strongly localized in any of the primary gaps, 
in full analogy with conventional periodic structures. 

We consider below wave transport and soliton motion 
in effectively discrete systems, such as arrays of optical 
waveguides or coupled-resonator structures in photonic 
crystals or arrays of traps for matter waves 20] where 
the wave transport is defined through tunneling between 
fundamental modes of the neighboring wells. This pro- 
cess can be described, in the tight-binding approxima- 
tion, by the discrete nonlinear Schrodinger (or Gross- 
Pitaevskii) equation for the normalized amplitudes E^, 

'i—^ \- Cn^n-l^n-l + C'n,n+l^n+l , 

dz (1) 

-^VnEn^T{n,\En?)En=^. 
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FIG. 1: Sketch of refractive index (or negative potential) in 
a lattice with a quasi-periodic modulation of (a) effective on- 
site potential, proportional to index contrast, and (b) coupUng 
between the lattice sites, inversely proportional to waveguide 
separation. 



Here n is the number of waveguide or potential well, 
z is the propagation distance or time, Cnm are cou- 
pling coefficients, Vn characterizes the detuning between 
the different sites, and defines the nonlinear response 
[^(n, 0) =0]. We note that the energy conservation law, 
P = l^nP = const leads to the requirement that the 
coupling coefficients are symmetric, Cn,n-i = Cn-i,n- 

The properties of discrete systems with a modulated 
potential Vn but homogeneous coupling C = const were 
extensively investigated with applications to the problem 
of electron transport and localization (see Refs. 0, 
|23,l2E07 and references therein). A special mathemat- 
ical theory was developed to demonstrate that the linear 
band-gap spectrum remains regular (absolutely contin- 
uous) for harmonic modulations of Vn below a certain 
threshold 25], whereas aperiodicity usually leads to the 
appearance of dynamically or critically localized states 
which are neither extended nor exponentially decaying 
inside the structure, and this inhibits efficient propa- 
gation of wave packets apart from specific resonances. 
There is one well-known example of aperiodic structures 
defined through modified Thue-Morse sequences, which 
supports extended modes and fundamental lattice soli- 
tons, however their motion is also accompanied by ra- 
diation due to mixed- frequency resonances and gap soli- 
tons were not identified [2^ . It remains an open problem 
whether there exist non-harmonic potentials which pro- 
duce regular wave spectrum. 

In engineered photonic structures such as waveguide 
arrays [3 there is a flexibility in designing both 
the on-site potential, and also the coupling coefficients. 
These can be controlled independently by adjusting the 
waveguide width and their separation, as schematically 
illustrated in Figs.^a,b). Recently demonstrated arrays 
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FIG. 2: Example of a lattice with quasi-periodic coupling 
and constant on-site potential [as sketched in Fig. db)] cal- 
culated with Eq. 0: (a) Coupling vs. lattice site; (b) Fourier 
spectrum of the coupling coefficient, shown with different ver- 
tical scale in top and bottom plots. 



of magnetic micro-traps 26] can be used to precisely de- 
sign, in a similar way, the trapping potentials for atomic 
condensates. We show in the following that specially de- 
signed modulation of inter-site coupling can provide the 
optimal conditions for linear and nonlinear wave trans- 
port in quasi-periodic systems. 

The model Eq. has an analog in the form of 
Ablowitz-Ladik system ^2} , 



dz 



+ (V^n-l + V^n+l)(C + iV^nP) = 0, 



(2) 



where = {i^n\ i^n\ - - -) is a vector function defin- 
ing the profiles of multiple interacting components. This 
is considered to be an integrable model [23, where soli- 
tons would interact elastically (without radiation losses) 
with other nonlinear waves, including extended quasi- 
periodic modes. Although nonlinear modes induce asym- 
metric coupling in Eqs. and therefore do not cor- 
respond to real structures, we show that by identify- 
ing the effective modulations of lattice coupling induced 
by quasi-periodic nonlinear waves in the framework of 
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Eq. we can then define equivalent 'optimal' modula- 
tion for the original Eq. (^). Then, we consider solutions 
of Eq. ^ in the form of stationary modulated waves 
i^n{z) = '0n(O) exp(zp^z), which profiles satisfy a set of 
difference equations, 

- Pni^n + + V^n+l)(C + |V^nP) = 0. (3) 

The waves with oscillating amplitudes tjjn induce, in gen- 
eral, a quasi-periodic lattjj3e with an effective modulation 
of coupling coefficients, Cn,n±i = + |V^nP- The lattice 
modulation can be controlled by choosing the number 
of modes, the values of pn, and the initial conditions i/jq 
and tjji . Dynamics of small-amplitude excitations in such 
lattices are governed by a linear equation with the mod- 
ulated coupling. 



-g^ + + ^n+l)C'n,n±l = 0. 



(4) 



Most importantly, Eqs. ^ are integrable l2q, and there- 
fore the band-gap spectrum of (j)n is completely regular, 
which is one of the key design requirements as outlined 
above. We now introduce a scaling transformation for 



the mode amplitude, (j)n ^ y Cn,n+i^ which preserves 
the regular band-gap spectrum and simultaneously trans- 
forms Eq. Q into a linearized Eq. with symmetric 
coupling found as 



(5) 



This is the key expression which defines an optimal cou- 
pling that can be implemented in real periodic structures. 

We show in Fig.[2fa) characteristic quasi-periodic mod- 
ulation of the inter-site coupling coefficient calculated us- 
ing Eqs. The Fourier spectrum of the coupling 
coefficient calculated over 1000 lattice sites contains 5 
dominant Fourier components [Fig.Efb, top)], which de- 
fine long-range quasi-periodic modulation. However, the 
bottom plot in Fig. Efb) clearly demonstrates that the 
spectrum structure is much more complicated, containing 
a continuous set of smaller components, which describe 
weaker aperiodic modulations of the lattice. 

In order to investigate the effect of deterministic aperi- 
odic modulations, we first analyze propagation of small- 
amplitude waves in the form En exp(i/3z), where the spec- 
tral parameter (3 has the meaning of the propagation con- 
stant of frequency for optical waves, or negative chemical 
potential for matter waves. For each value of (3, linearized 
Eq. admits two linearly-independent solutions, 
and E~ . Since the linear spectrum is regular being con- 
nected to integrable equations, are either both ex- 
tended indicating the presence of transmission band, or 
exponentially localized when /3 is inside a band-gap. In 
the latter case, we can choose the mode indices such that 
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FIG. 3: Top: Band-gap diagram shown as inverse localiza- 
tion length vs. the propagation constant. Bottom: Moving 
discrete and gap solitons corresponding to the total internal 
reflection and Bragg-reflection gaps as indicated by arrow. 
Solitons are excited by a single tilted beam (discrete) or two 
tilted beams at the Bragg angle. The lattice parameters cor- 
respond to Fig. 121 
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FIG. 4: Band-gap structure and soliton dynamics in a pure 
"harmonic" lattice, which coupling spectrum contains only 
the 5 dominant peaks shown in Fig.|2Ib, top), without aperi- 
odic contributions. Input conditions are the same as in Fig.|3 



lim^^ioo = 0? can be demonstrated using the 

transfer-matrix approach that their exponential de- 
cay rate is the same. Numerically, we calculate the in- 
verse localization length as follows. 



Hi{/3) = lim log 



max \En\ 

n/2,n 



max I En I 
l,n/2 
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This dependence defines the hnear spectrum: ^ 
for extended waves inside the bands, and k > cor- 
respond to band-gaps, see an example in Fig. [Sftop). 
There is a fundamental total-internal-reflection gap, as 
well as additional gaps corresponding to Bragg-reflection 
resonances at the dominant modulation frequencies [cf. 
Fig. [Hb,top)]. Each of these gaps can support bright 
solitons, that may move through the quasi-periodic lat- 
tice in a full analogy with periodic structures 0, 
El El E 113, see examples in Figs. Ofbottom). These 
simulations are performed for a self-focusing Kerr-type 
nonlinearity, = |£^np7 however similar results are ex- 
pected for other nonlinear dependencies. We note that 
some radiation is observed at the initial stage due to non- 
ideal input conditions modelling one-beam 14] and two- 
beam 0, experimental arrangements, however radi- 
ation losses of propagating solitons are practically negli- 
gible, resembling radiation in periodic lattices j^fll that 
remains exponentially small in a broad parameter region. 

For comparison, we also consider a 'pure' quasi- 
periodic lattice, which modulation spectrum contains 
only the five dominant components shown in Fig. EJb, 
top). Its band-gap spectrum [Fig. ^top)] contains the 
same principal gaps as observed in Fig. inijtop). How- 
ever, there also appear multiple additional mini-gaps due 
to mixed-frequency resonances, which are no longer sup- 
pressed by aperiodic modulations. We perform numer- 
ical simulations using the same input conditions as in 
Fig. Efbottom), and observe that presence of mini-gaps 
strongly affects the soliton motion, see Fig. Etbottom). 
Indeed, the motion of discrete soliton is arrested as it 
becomes trapped, similar to effects in disordered sys- 
tems 5 the gap solitons emit strong radiation dur- 
ing their propagation. It is most remarkable that small 
aperiodic modulations have such a dramatic effect on lin- 
ear wave propagation and soliton dynamics. 

In conclusion, we have demonstrated that linear wave 
propagation and soliton motion in lattices with exactly 
quasi-periodic modulation of their parameters can be af- 
fected by Bragg resonances at mixed modulation frequen- 
cies, leading to soliton trapping or strong radiation losses. 
However, these negative effects can be suppressed in lat- 
tices with additional aperiodic modulations, and we have 
presented a systematic analytical method for calculating 
the parameters of such optimized lattices. These results 
can be used, for example, to design photonic structures, 
such as waveguide arrays, and arrays of traps for atomic 
condensates. 
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